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The Cosmic Defect theory has been confronted with four observational constraints: primordial nuclear 
species abundances emerging from the big bang nucleosynthesis; large scale structure formation in the 
universe; cosmic microwave background acoustic scale; luminosity distances of type la supemovae. 
The test has been based on a statistical analysis of the a posteriori probabilities for three parameters 
of the theory. The result has been quite satisfactory and such that the peiformance of the theory is not 
distinguishable from the one of the ACDM theory. The use of the optimal values of the parameters for 
the calculation of the Hubble constant and the age of the universe confirms the compatibility of the 
Cosmic Defect approach with observations. 
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1. Introduction 

From the very beginning of the theory of relativity both in its special form (SR) and in its 
general form (GR) a problem of interpretation of the nature of space-time has been present. 
This problem, in most cases, is implicit rather than explicit, but it is there. Is space-time a 
sort of field? which would be perfect for people trying to quantize the gravitational inter- 
action. This cannot be, however, since, in general, fields are described in space-time; what 



1 
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would then be the background of space-time? What in any case is granted is that space-time 
is indeed something physical, not a mere mathematical artifact. Space-time interacts with 
matter as the Einstein equations tell us, so even not being matter it has properties on its 
own. 

On another side we know that space-time is a curved four-dimensional manifold en- 
dowed with Lorentzian signature, and we also know that space, at the cosmic scale, ap- 
pears to be expanding with a typical symmetry: the Robertson-Walker (RW) s ymm etry. 
Even more: some pieces of observation tell us that the expansion is accelerated 'JJ^I 'YYas, 
acceleration is not evidently due to matter, so commonly it is ascribed to "something else" 
which is du bbed in various ways acc ording to different theories, but is mostly known as 
dark energy El 4 | 5 1 6 | 7 | 8 | 9 | 10 | 1 1 | 12 | 13 | w^jj^j^ jg dscnk energy? There are as many answers as 
there are theories, but it is difficult to say that the situation is clear and that each interpreta- 
tion of dark energy is physically well motivated. 

We refer here to the Cosmic Defect (CD) theory. It describes space-time on the basis 
of the existing analogy with usual three-dimensional elastically deformable continua; the 
classical theory is of course extended to four dimensions and endowed with the appropri- 
ate signature. All features of GR are preserved. The theory is exposed in Eland will be 
shortly reviewed in the next section. There are two relevant features with CD: a) the global 
symmetry of the universe is due to a defect in the texture of space-time working as defects 
do in material continua; b) the deformation of the manifold, due both to defects and to the 
presence of matter is expressed by a strain tensor which coincides with the non-trivial part 
of the metric tensor To the strain tensor an elastic potential energy density corresponds, 
whose presence affects the global behavior of space-time and is responsible for the acceler- 
ated expansion. The theory is in our view appealing, but of course it must be tested on facts. 
We already used it to work out the luminosity/distance curve of type la supemovae with 
good results however the fit of the luminosity of Snla's is not difficult to be obtained 
by a variety of different theories so it is not really discriminating among them; this is the 
reason why we decided to test CD against more constraints. One such constraint concerns 
the cosmological nucleosynthesis (Big Bang Nucleosynthesis: BBN): the expansion law 
of the universe must be such that the appropriate pressure and temperature conditions last 
long enough (and not more) to produce He and D in the percentages we observe or deduce 
from indirect observations K 



We considered two other main tests: the large scale structure (LSS) formation, which 
probes the universe history at the time of matter-radiation equality, and the acoustic scale 
of the cosmic microwave background (CMB), which gives an independent constraint at 
the recombination era. This allowed us to compare the performance of the CD theory with 
other theories and especially with ACDM. After exposing the method that we have adopted 
in order to validate the constraints, we shall verify also the acceptability of the estimates 
obtained from CD for the Hubble constant and the age of the universe. We shall see once 
more that the results are more than acceptable and we shall discuss them in the concluding 



^More light elements would enter a more refined and detailed analysis but the two considered in the text are 
enough for our puiposes. 
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section adding some consideration regarding the Solar system environment. 
2. Review of the CD theory 

The core idea of the Cosmic Defect theory is that the actual space-time manifold with its 
global and local curvature behaves as a four-dimensional elastic continuum, so that one may 
think that the natural situation is obtained introducing strain in an initially flat (Minkowski) 
manifold. The properties of the strained manifold are expressed in terms of two parameters, 
which are the Lame coefficients of space-time, A and /i. The details of the theory may be 
found in 231 jigre we only recollect that the global RW symmetry is assumed to be induced 
by a cosmic defect and that the actual metric tensor is composed of two contributions: 

gtif = Vfj.i' + 2e,ii/, (1) 

where ij^i, is the Minkowski metric tensor and e^^ is the strain tensor According to the 
analogy with a deformed elastic material we expect the strain to be associated with an 
elastic deformation energy written as: 

We^^Xe^ + ^^e^pe"^. (2) 

Now e — e^is the trace of the strain tensor The energy in (|2) allows to write a new action 
integral including both space-time and matter/energy in the form 




^matter IS the Lagraugiau density of ordinary matter/energy and the rest has the usual mean- 
ing. The "elastic" potential energy term belongs to geometry, i.e. space-time, even though 
it looks like some matter contribution. Considered from the field theoretical viewpoint this 
new contribution implies a "mass" associated with the gravitational interaction; usually this 
is said as the graviton being massive, which fact has relevant consequences when studying 
the propagation of gravity and, in particular, of gravitational waves In any 
case our conceptual framework is entirely classical so that, p roperly speaking, there are no 
gravitons, but rather it turns out that gravity has a finite ranget^SIE! We are not discussing 
the issue in this paper; we only remark that the numerical values we find after our tests and 
fits to observation tell us that the effects of the space-time strain are totally negligible at the 
scale of the Solar system. 

As we have seen, everything depends on the strain tensor, which in turn depends on 
the way events on the reference manifold are associated to their corresponding events in 
the natural manifold. Actually there are infinite possible ways to get a given final situation 
starting from a flat initial one. This apparent freedom of choice has indeed a physical mean- 
ing since our manifolds are physical. So different choices correspond to different strains and 
the Hamilton principle permits to identify the least strain configuration. In practice what 
can be seen as a gauge freedom appears as a gauge function in the line element of the 
flat reference written using the coordinates of the curved natural manifold. We expect the 
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different gauge choices not to affect the general laws, such as the expansion law in the cos- 
mological case, rather to point out different properties of space-time, i.e. different values 
of its Lame coefficients. 

Applying all the above to the RW symmetry in (|3]l one is left with two free functions. 
One of them is physical: the scale factor a(r), depending on the cosmic time r. The other is 
the mentioned gauge freedom expressing the difference between cosmic time and the time 
coordinate on the flat reference manifold; in other words, using the cosmic coordinates, the 
flat line element of the reference manifold is written 



ds^ = b\T)dT^ - dx' - dy' - dz^ (4) 
and hir) is the mentioned gauge function. 

From the RW Une element and from (|4]l it is straightforward to write the non-zero 
elements of the strain tensor, which are 

Coo - 



2 

l-a2 



(5) 

Introducing the explicit form of the strain tensor into the Lagrangian of our problem, 
then applying the usual variational principle, one arrives to a pair of equations: one is solved 
to fix the gauge i.e. to give the explicit dependence of h on a; the other, after one integration, 
gives the expansion law of the universe, which is implicitly written in the formula for the 
Hubble parameterEl 

^ ^ a _ /SLA (1 
a 




8k 



1/2 



■ 3_g-(l + 2)'[Pm0 + Pr0(l + z)]| . (6) 

Dots are derivatives with respect to cosmic time; z is the redshift factor; ao is the value 
of the present cosmic scale factor; k ~ IGttG/c^ is the coupling constant between mat- 
ter/energy and geometry; pmo is the present mass density of the universe and pio is the 
present radiation energy density. The B parameter is the bulk modulus of space-time and 
its dependence on the Lame coefficients is obtained from the gauge condition. It is: 

that works for any combination of the Lame coefficients excepting A = — 2/i0 

In the approach corresponding to (|6) the content of the universe is simply dust and 
radiation. Formula (|6) is the starting point for the elaborations in the rest of the paper 



''The expression found inJ^has been obtained forcing the gauge to be 6 = 1 which is only consistent with the 
condition A = 0. 
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3. Observational constraints 

Observational cosmology allows us to put strong constraints on the expansion history of the 
universe over a very large span of time. Here, we want to test if the scale factor evolution in 
(|6]l is compatible with such history. In order to perform a conservative analysis, we looked 
at geometrical tests based on a minimum set of hypotheses. We did not consider evidence 
more heavily dependent on the fraction of baryonic matter (such as peaks in the galaxy 
power spectrum) or on the present linear-theory mass dispersion as- Constraints from the 
cosmic microwave background (CMB), which require modeling of additional parameters 
for spatial curvature, amplitude and slope for the tensor spectrum and optical depth to last 
scattering, were not considered too. We have however estimated the CMB acoustic scale 
which can show whether our theory behaves differently from ACDM and other theories or 
not. 

As standard relativistic contributions, we considered a photon background at T = 
2.728 kI^SI together with three species of massless neutrinos. 

3.1. Nucleosynthesis 

At very early time, the matter/energy budget is radiation dominated and the Hubble factor 
can be approximated as 



The strain energy is not supposed to affect the cross-sections of local nuclear interactions, 
whereas, as far as the expansion rate is considered, it effectively boosts the radiation density 
by a factor ATeoost = (1 + B/Bag). Then, the strain of space-time enters nucleosynthesis 
processes only through a modified expansion. Due to the increase in the overall rate of 
expansion, neutron freeze-out happens earlier, raising the final helium abundance. Non- 
standard boost factors (A'soost 7^ 1) can then be tightly constrained by fitting measured 
abundances. 

Several methods have been proposed to determine the ''iJe mass fraction, Y^l^J-lwith 
systematic errors being the main responsible for the spread determination. To account for 
this, EH performed an analysis of several values in literature, ending up with an estimate of 
Yp — 0.250 ± 0.003. This range of values is still fully compatible with other more recent 
estimates ^221 ^ve still refer to the conservative analysis in '211 The related constraint on 
the boost factor is XBoost = 1-025 ± 0.015l^, from which we can write a constraint as 




(8) 



with 



(9) 




(10) 
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3.2. Structure formation 

Matter perturbations can not grow in a universe expanding as a radiation-dominated back- 
ground. The space-time strain effectively increases the radiation-like density and matter 
dominance is postponed. The effective boost Xsoost affects the scale of the particle hori- 
zon at the equality epoch, z^q ~ SlSoES On the other hand, the Newtonian limit of gravity 
is still fulfilled in presence of defects El so that, in a CD model, perturbation growth is 
affected mainly through the modified expansion rate of the background. The horizon at 
the equality is imprinted in the matter transfer function. The constraint from large scale 
structure (LSS) becomes 

(^mo/l) apparent = ^'^Boost (^mo/l)t,.uo " 1) 

where h is the Hubble constant Hq in units of 100 km s^^Mpc"^ and fimo is the matter 
density in units of the critical density pcr = SHq / {8ttG). 

We considered the final analysis from the 2dF Galaxy Redshift SurveyE3which found 
(^2i„o/i)apparcnt = 0.168± 0.016. The fitting procedure inESIwas performed trying to Umit 
the number of parameters so that you can reliably use their result combined with Eq. (fTTT i 
to get constraints on extra-physics. We share the same assumption on the index of the 
primordial power spectrum (n — 1). The related constraint on the cosmological parameters 
of the CD theory can then be written as 

, _ / ^Boostano/^- 0-168 y 

^LSS - ^ j ■ (12) 

3.3. Cosmic Microwave Background 

The temperature power spectrum of CMB is sensitive to the physics at the decoupling 
epoch, zls ~ 1090, as well as the expansion history between now and the last scattering 
surface. Some of the main features of the spectrum can b e summ arized by a number of 
parameters. Among them, we consider the acoustic scale 

?A = (l + .Lsh^^^, (13) 
fs(2Ls) 

where rg is the sound horizon at recombination and Da is the angular diameter distance 
to the last scattering surface. The cosmological model we are testing affects the power 
spectrum in two main ways. Firstly, the radiation boost Xeoost affects the location of the 
acoustic peaks, which depends on the expansion factor at the matter-radiation equality, Ueq- 
The formula for the sound horizon is then the same as for the ACDM model ESEZl with 
the only caveat that now — ^BoostPro/PmO- Secondly, the angular diameter distance 
depends on the expansion history, 

(1 + ZLs) 7o H{z) 

The values determined for I a are quite model-independent, so that we can consider the 
latest results from WMAP-7 ^3. We take l^*"" = 302.69 ± 0.76 ± 1.00, where the first 



12,2013 21:47 WSPC/INSTRUCTION FILE defects ijmpd vS 



Constraining CD theory 7 



error is the statistical error in a ACDM model and the second error gives an estimate of 
the uncertainty connected to the model . Since these errors are independent they can be 
added in quadrature. The statistical constraint then reads 

2 f UjB, n^o, ggj- 302.69 V 
XcMB^[ 1:26 ) ■ ^^^^ 



3.4. Supernovae 



Supernovae (SNe) of type la are supposed to be the best tracers of the recent expansion 
19 01 

history of the universe. 1^ provided a sample of 307 SNe for cosmological analysis. As 
usual, we can compare the measured to the predicted distance moduli. 



d{z) = 25 + 5 log 



10 



(1 + ^) 



(c/Mpc) ; 



The related term is 



XsNc 



E 



d{zi) 



Ad,: 



(16) 



(17) 



where d, is the measured distance modulus at 



4. Statistical analysis 




Fig. 1. Posterior probability density functions for the parameters pmO (left panel), Bag (middle) and B (right). 
Units are as in Table[T] 



Table 1. Estimated parameter values. The maximum likelihood estimates are 
reported in brackets. 

B (10--'^^m-^) I pmo(10-29gcm-3) | B-g\l0^^m'^) 
2.22 (2.22) ±0.06 I 0.260 (0.258) ± 0.009 I 0.011 (0.009) ± 0.006 



We performed the statistical analysis with standard Bayesian 

methods™ 

The Bayes 

theorem states that 

p(P|d)(x£(P|d)p(P), (18) 
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where p(P|d) is the posterior probability of the parameters P given the data d, £(P|d) 
is the likelihood of the data given the model parameters and p(P) is the prior probability 
distribution for the model parameters. 

For our analysis, three model parameters are enough to describe the CD theory. The 
discussed constraints do not require us to distinguish between baryonic and dark matter, so 
that in the following analysis we can consider a single parameter density for the dust-like 
matter, pmo. Cosmic defects characterize the B parameter. Finally, the present value of the 
scale factor is described in terms of Bao ■ For convenience, we actually used as parameter 
the inverse of Ba,,- Finally, P = {pmo, B, B^^}. 

In order to investigate the parameter space for the CD model, we have to write the 
likelihood. The constraints discussed in Section|3]allow us to write the total as 

X^{PnM,B,B^^\d) = xInc + xIbN +XLSS +XCMB- (19) 

The related likelihood is then C cc e^^ 

As a priori information, we considered flat distributions for each parameter The rela- 
tivistic energy density was fixed at p^o ~ 7.8 x 10^'^'' g/cm . 

The parameter space was explored with Monte Carlo Markov chains methods. We run 
four chains, each one with lO** samples. Convergence criteria are safely satisfied, with the 
Gelman and Rubin ratio ^1 29 jjgjjjg < i .003 for each parameter Results are summarized 
in Fig. [Hand Tabled 

The posterior probability densities (pdfs) are plotted in Fig.[Tl Distribution s have been 
smoothed with a Gaussian adaptive kernel with reflective boundary conditions EZEl] -pjjg 
pdfs are very regular and single peaked. Parameter estimates, computed as mean and stan- 
dard deviation of the final marginalized distributions, are reported in Table [T] We also per- 
formed a standard maximum likelihood analysis. The maximum was found with a downhill 
simplex algorithm . As expected for single-peaked distributions, Bayesian estimates co- 
incide almost perfectly with results from the maximum likelihood analysis. 

The elastic parameter B is mainly determined by late time evolution, when the space- 
time strain energy propels acceleration, whereas Bag is pretty much constrained by the 
BBN, the constraint from LSS being weaker. The matter density is well determined thanks 
to LSS and SNe. The density parameter reads fl^o = 0.28 ± 0.01. 

Today's scale parameter is not well constrained by the data. LSS and BBN require qq 
to be large enough so that the boost factor Xboost is small at early times. We got gq > 19.9 
at the 99% confidence level (CL) and Xboost = 1-024 ± 0.013 with Xboost < 1-053 at 
the 99% CL. On the other hand, SNe data are compatible also with much smaller values of 
order unity for oq. Note that the tail at large values for the a posteriori pdf of oq is strongly 
dependent on the chosen prior on B^^, whereas the lower limit is not. 

5. Compatibility checks 

5.1. Some observational predictions 

To further test the performance of the cosmic defect theory, we can test some predic- 
tions based on the above statistical analysis. The inferred Hubble parameter turns out to 
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be Hq — 70.2 ± 0.5 km s^^Mpc^^, in very good agreement with the local estimate of 
73 ± 2 ± 4 km s^^Mpc^^ based on high precision distance determination methods, SSI 
The estimated age of the universe is to = 13. 7 ± 0.1 Gy in good agreement with lower 
limits obtained by the ages of the oldest globular clusters (12.6^2 2 95% CL*^^ and 
radioactive dating (12.5 ± 3 GyES). 

5.2. Comparison to ACDM 

Table 2. and BIC values for the CD model versus ACDM. 



BIC 

Constraints CD ACDM CD ACDM 



SNe 


310.3 


311.9 


327.5 


323.4 


SNe+LSS 


311.9 


313.8 


329.1 


325.3 


SNe+LSS+CMB 


314.1 


314.1 


331.3 


325.6 


SNe+LSS+BBN 


312.3 




330.8 




SNe+LSS+BBN+CMB 


315.0 




332.3 





CD model predictions for ilmo and Hq are in good agreement with results for a 
flat ACDM model Let us now discuss how these two competing models compare 
under a statistical point of view. The quickest way to model comparison is to exploit 
some information criterium which takes into account how well a model fits the data 
versus its complexity ED. in particular, the Bayesian information criterium is defined as 
BIC — + Npar log Nconst, whcrc is the total for the model, Npar is the num- 
ber of parameters and Nconst is the number of observational constraints. The best model 
minimizes the BIC. 

The CD model has three free parameters; two parameters, {Hq, ilmo}, characterize 
the flat ACDM model, performances are listed in Table |2] Comparison was performed 
under hypotheses favorable to ACDM. We did not use the BBN constraint, since the stud- 
ied ACDM model lacks of any additional relativistic species to fit the abundances. This 
circumstance further penalizes the model with the larger number of parameters. However, 

and BIC values are too close to prefer one model over another 

A full treatment would require the computation of Bayesian evidences. However, due 
to the very similar likelihood functions, the estimated evidences would be highly dependent 
on the assumed priors, which is not really informative. 

6. Conclusions 

We have submitted the CD theory to a consistency test with respect to three main cosmolog- 
ical constraints: primordial nuclear isotopic abundances; large scale structure formation in 
the universe; luminosity distances of type la supemovae. One relevant quantity, the sound 
horizon, deducible from the CMB anisotropy data has also been evaluated. According to 
the CD theory space-time is endowed with a strain energy density, whose presence in the 
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Lagrangian density of space-time affects the expansion law of the universe and produces 
the present accelerated expansion. In early times the effect of the strain term looks like 
the one of radiation, however it shows up only in the expansion rate without directly af- 
fecting either the interaction between radiation (photons and neutrinos) and particles or the 
cross sections of the nuclear reactions. In practice the relevance of the parameters of the 
CD theory appears in a boost factor. Similar considerations hold for the large scale struc- 
tures formation as well as for the sound horizon of the CMB at last scattering. Later, the 
accelerated expansion, which shows up as a dimming of the Snla's, is essentially deter- 
mined by the bulk modulus of space-time treated as an elastic four-dimensional manifold 
with Lorentzian signature. The constraints posed by observation may be translated into 
constraints for the values of three parameters of the CD theory. On these bases a statistical 
analysis has been done on the a posteriori probabilities of the values of the parameters. By 
these means we have found: a) that CD is compatible with the observational constraints; b) 
that the results are as good as the ones obtained from the ACDM theory. These conclusions 
are strengthened by the use of the optimal values of the parameters for determining the 
Hubble constant and the age of the universe: the values we obtain are in good agreement 
with the currently accepted ones. 

The tests we performed on CMB and LSS implicitly assumed a Newtonian approach 
for the density perturbation growth. We can check a posteriori that such assumption holds. 
The LSS test in Sect. [3] was based on the main bend visible in the matter transfer function, 
which arises because the universe expansion is radiation dominated at early times whereas 
fluctuations in the matter can only grow if dark matter and radiation fall together Similar 
considerations stand for the position of the acoustic peak in the CMB spectrum. Both tests 
then exploit a distinctive geometric feature in the respective spectra. 

As for the expansion law of the universe, we used the exact result of the CD theory. For 
the gravitational potential, we assumed a locally Newtonian potential, i.e., a solution of the 
standard Poissonian equation. This is fully justified given the size of the Lame coefficients 
of space-time, which determine the relevant scale for the CD theory. 

The test with the SNe is independent of the growth of perturbations and relies only 
on the expansion history. It constrains B, whose size in turn gives that of A and ji. The 
SNe analysis unambigously shows that the order of magnitude of the Lame coefficients of 
space-time must be ^ lO^^^m^^, see also^l. On a dimensional analysis, local deviation 
from the Newtonian potential near a mass m must be of the kind of Ar^, fjLr^ or Amr, fjLmr 
or higher order corrections. Even on the present scale of a galaxy cluster (to ^ 10^'* Mq, 
r ^ 1 Mpc), the Newtonian term m/r is a couple of orders of magnitude larger than the 
CD correction. It is then safe to consider a standard Newtonian potential in the test we 
performed. The basic version of our tests on CMB and LSS does not then require a full 
treatment of perturbation growth in the CD theory. 

Even if the CD theory provides a novel theoretical framework, the above considerations 
are the same generally done for the cosmological constant. Despite the very different origin, 
A is estimated to be of the order of ^ lO^^^m^^ too. Its contribution to the local potential 
is neglected in first approximation when considering density perturbations and the simplest 
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One more remark we can make concerns the possible Solar System effects of the strain 
of space-time. A full discussion of this issue requires the explicit solution of the Einstein 
equations with strain in a spherically symmetric stationary case, which is beyond the scope 
of the present paper Let us again compare B times the square of the distance from the 
center of the system to the equally dimensionless quantity m/r that tells us how strong are 
the gravitational effects for instance in the Solar system. In the environment of any stellar 
system and for distances of a few astronomical units m/r is of the order of ~ 10~^°~10~^. 
Considering the result in Table[T]we see that, in the same distance range, it is Br^ ^ 10^^*, 
i.e. some 20 orders of magnitude below the Newtonian mass terms. For this reason we 
expect the effects of the strain of space-time to show up at cosmic scales or at least at very 
large scales, but not within the Solar system. 

Two appealing features of our approach stand out. On the theoretical side, CD theory 
is a paradigm based on a physical interpretation, at least by analogy, of the behavior and 
properties of space-time and not simply a mathematically consistent frame where facts fit 
without a corresponding intuitive picture of the situation. On the observational side, the 
viability of CD is as good as the one of ACDM. The positive tests described in this paper 
steer us to continue to work out all consequences implied in the Cosmic Defect theory, in 
quest not only of an alternative picture, but also of some peculiar property allowing really 
to discriminate one paradigm from the others. 
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